Long range Josephson coupling through ferromagnetic graphene 
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We study the Josephson effect in graphene-based ballistic superconductor-ferromagnet- 
superconductor (SFS) junctions. We find an oscillatory Josephson coupling I c Rn of F graphene 
whose amplitude is nonvanishing for a half-metallic graphene, increases for the exchange fields h 
above the Fermi energy Ef and shows only a slow damping at strong exchange fields h S> Ef- We 
interpret this long range Josephson coupling as the result of the exchange mediated Andreev-Klein 
process at FS interfaces which enhances the induced antiparallel-spin superconducting correlations 
in F graphene by increasing h above Ef- We further demonstrate the existence of regular tem- 
perature induced transitions between and n couplings in the plane of T and h where the phase 
boundaries have distinct shapes at the two regimes of h below and above Ef- 
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Graphene, the two dimensional (2D) solid of carbon 
atoms with honeycomb lattice structure, shows unique 
properties due to its peculiar gapless semiconducting 
band structurei^^. The conduction and valence bands 
in graphene have the conical form at low energies with 
the apexes of the cones touching each other at the cor- 
ners of the hexagonal first Brillouin zone which deter- 
mine two non-equivalent valleys in the band structure. 
The charge carrier type (electron- like (n) or hole-like (p)) 
and its density can be tuned by means of electrical gates 
or doping of underlying substrate. An important aspect 
in graphene is the connection between its specific band 
structure and the pseudo-spin which characterizes the 
relative amplitude of electron wave function in two dis- 
tinct trigonal sublattices of the hexagonal structure. This 
has caused that the charge carriers in graphene behave 
as 2D massless Dirac fermions with a pseudo-relativistic 
chiral property*. Already anomality of variety of phe- 
nomena including quantum Hall effect, Andreev reflec- 
tion (AR)»£, and Josephson effeclj^ in graphene have 
been demonstrated. Here we report on the peculiar- 
ity of Josephson effect in a graphene superconductor- 
ferromagnet-superconductor (SFS) junction, which arises 
from such a Dirac-like spectrum with chirality. We find 
that a weakly doped graphene F contact can support a 
long ranged opposite-spin supercurrent which persists at 
strong spin-splitting exchange fields, in striking contrast 
to the behavior of the Josephson current in common SFS 
junctions^. 

Superconducting correlations can propagate through a 
mesoscopic normal metal (N) contact between two super- 
conductors (S) via the process of AR at the NS-interfaces 
in which the subgap electron and hole excitations with 
opposite spin directions are converted to each other—. 
Successive AR at the two NS interfaces and the coher- 
ent propagation of the excitations between these reflec- 
tions leads to the formation of the so called Andreev 
bound states which can carry a supercurrent. The result- 
ing Josephson effect, characterized by the critical (maxi- 
mum) supercurrent I c and a relation with the phase dif- 
ference tp between superconducting order parameters of 
the two superconductors, is well established in a vari- 



ety of SNS structures^. In an SFS junction, due to 
the exchange correlations field h, a momentum change 
of 2h/vF between Andreev correlated electron-hole is in- 
duced which results in a damped oscillatory variation of 
I c with the length of F-contact L. As the result of the 
J c -oscillations an SFS structure can transform into the 
so-called 7r-junction in which the ground-state phase dif- 
ference between two superconductors is tt instead of 02ii£. 
The damping of I c occurs over the magnetic coherence 
length £/j which is ~ hvp/h for a ballistic F 13 i 14 . This 
makes the Josephson coupling in F junctions rather short 
ranged as compared to SNS systems in which the Joseph- 
son coupling persists over much longer lengths of order of 
the normal-metal coherence length £jv = hvp /kuT 9 (nor- 
mally h is much larger than the superconducting gap A). 
In particular for a half-metal F contact with h > Ep 
there is no Josephson coupling for a sizable contact of 
length L > Xp^. 

In this work, we demonstrate unusual features of the 
exchange-induced / c -oscillations and the corresponding 
— it transitions in a ballistic F-graphene Josephson con- 
tact between two highly doped superconducting regions 
(see Fig. [1]). We show that while in the regime of h < Ep 
the amplitude of the critical current shows a monotonic 
damping with the exchange field, for the higher exchange 
fields h > Ep it develops drastically different behavior. 
For a half- metal F with h — Ep, we find that the Joseph- 
son coupling I c Rn (Rn being the normal state resistance 
of the junction) has a non-vanishing value, in spite of the 
vanishing density of states for spin-down electrons. In- 
terestingly, this finite Josephson coupling is resulted from 
particular Andreev bound states in which spin-up propa- 
gating excitations and spin-down evanescent excitations 
are involved. These mixed evanescent-propagating states 
can have significant contribution in the supercurrent due 
to the chiral nature of the carriers in F graphen o 7 ! 15 . 

More surprisingly, we find that for exchange fields 
above the Fermi energy h > Ep the coupling I c Rn in- 
creases above its half-metal value and shows damping 
only at strong exchange fields h 3> Ep with a rate which 
is much lower than that of the regime of h < Ep. We 
explain this long range Josephson effect in terms of su- 
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perconducting correlations between a n-type excitation 
from the spin-up conduction subband and a p-type exci- 
tation from the spin-down valence subband in F (see Fig. 
[TJd). For h > Ep these two types of excitations are cou- 
pled at the FS-interfaces via a peculiar Andreev process 
which is accompanied by a Klein tunneling through the 
exchange field p — n barrier 16 . It has been found that this 
spin Andreev-Klein process leads to an enhancement of 
the amplitude of AR and the resulting subgap conduc- 
tance of FS junctions with the exchange field. In the SFS 
structure the corresponding Andreev-Klein bound states 
are responsible for the long range proximity effect. We 
further demonstrate the existence of the temperature- 
induced regular — tt transitions by presenting phase 
diagram in T/T c and h/ Ef plane where the boundaries 
of — 7r phases have different forms in two regimes of 
h < Ep and h > Ep. 

To be specific, we consider a ballistic F-graphene strip 
of length L smaller than the superconducting coherence 
length £ = Hvf/A which connects two S electrodes (see 
Fig. []}. Highly doped superconducting regions can be 
produced by depositing superconducting metallic elec- 
trodes on top of the graphene sheet 8 . A graphene SFS 
structure similar to our setup has been studied by Linder 
et al<^ considering certain values for the Fermi energy in 
F and in the electrodes. They concentrated on the ex- 
istence of a large residual supercurrent in the points of 
— 7r transitions at T = 0. Here we consider a more 
realistic model of highly doped S electrodes and cover 
full range of the key parameter h/ Ep to find the above 
mentioned long range Josephson coupling whose under- 
lying mechanism will be explained in the following. We 
take the Fermi wavelength A^ s of S-electrodes to be very 
smaller than the superconducting coherence length £ and 
the Fermi wavelengths of the two spin subbands in F- 
graphene Afv- By the first condition mean field theory 
of superconductivity will be justified and by the second 
we can neglect the spatial variation of the superconduct- 
ing order parameter A(x) in the superconductors close 
to the FS interfaces. Thus A(x) has the constant val- 
ues A exp(±i</j/2) in the left and right superconductors, 
respectively, and vanishes identically in F. 

The superconducting correlations between a spin a 
electron excitation of wave function u a and the spin a 
hole excitation of wave function v s can be described by 
Dirac-Bogoliubov-de Gennes (DBdG)£ equation which, 
in the presence of an exchange field, reads: 

{ 6 '-£*- c Er -i -.)(::)-• « 

Here H a — —ihvp(d x a x + d y a y ) — ah and A = Aoo 
are the spin-er single-electron Dirac Hamiltonian and the 
superconducting pair potential, respectively and e is the 
excitation energy. The wave functions u a and v s are 
two-component spinors of the form (V'i, ^2) and <Xj (i = 
0, x, y, z) are Pauli matrices, all operating in the space of 
two sublattices (pseudo-spin) of the honeycomb lattice. 




FIG. 1: (Color online) (a-b) Schematic of the graphene SFS 
junction and the configuration (being electron-like (n) and 
hole-like (p)) of the spin- up and spin-down subbands for two 
regimes of (a) h < Ef and (b) h > Ef- The orientation of 
the wave vectors and velocity vectors of Andreev correlated 
electron-holes (for a normal incidence of electron to the in- 
terface) is also shown. For h < Ef the retro reflected hole 
has antiparallel momentum and velocity, whereas for h > Ef 
Andreev reflection is specular with the hole having parallel 
momentum and velocity, (c) Phase diagram of the — n tran- 
sition of the junction of EfL/Hvf = 10 around h = Ef- The 
boundaries between and n phases for h > E and h < E are 
the mirror form of each other. 

Inside F the solutions of DBdG equation JT]) are elec- 
tron and hole-like wave functions which are classified by a 
2D wave vector k CT = (fc CT , q) with the energy- momentum 
relation e a = ftUi?|k CT |. For a finite width W the trans- 
verse momentum is quantized {q n = (n + l/2)n/W) by 
imposing the infinite mass boundary conditions at the 
edgesi^. At the Fermi level e = for a spin direction 
a and a given q n there are two electron and two hole 
states with the wave functions u„ = v a = exp(±ik a x + 
iqy)(l, ± exp(±icv CT )) which are characterized by longitu- 
dinal momentum k a = \J fclT — q 2 and the propagation 
angle a a = arcsin(g/ kp a ) (kpu = (Ep + ah) / '(fivp) being 
the Fermi wave vector of spin a subband). 

The solutions of Eq. |T]) inside S (h = 0) are rather 
mixed electron-hole excitations, the so called Dirac- 
Bogoliubov quasiparticles. Assuming ideal FS contacts 
the electron-hole conversion can be described by a bound- 
ary condition between electron and hole wave functions 
which for the left and right interfaces, respectively, has 
the formsi, 

u a = e w/2 +l /3„. CT y _ i p = arccos ( e / A ), (2) 

where n is the unit vector perpendicular to a FS interface 
pointing from F to S. 

Introducing the normal-state transmission coefficient 
of spin a quasiparticles through the junction as t a — 
\t a \ exp(ir) a ) — (cosj a — zsin7 (T /cosQ; (T ) _1 with j a — 
k a L and imposing the conditions {2| at the two FS 
boundaries (x — 0, L), we obtain the following result for 
the energy of the spin a Andreev bound state 

e ff = Acos[(0(0)+7te-Tfe)/2], (3) 
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FIG. 2: (Color online) Dependence of zero temperature 
Josephson coupling I c Rn on the exchange energy h/ Ef scaled 
with Fermi energy for the doping EfL/Hvf = 10. Cusp-like 
variations indicate — tt transitions with a period of order 
hL/hvF- The coupling has a nonzero value for Ii/Ef = 1, 
develops a smooth maximum for K/Ef > 1 and shows a slow 
damping at strong exchange fields. For }i/Ef > 1 the scale of 
H/Ef is logarithmic to clarify the slow decrease in the Joseph- 
son coupling. Inset compares the same dependence at T = 
and T = 0.99T C . 



where cos 6 — \t a t s \ [cos 4> + tana ff tan a s sin7 CT sin7 ff ]. 
For a short junction of L <C £ only the Andreev bound 
states with energies |e| < Ao have the main contribution 
to the supercurrent. At temperature T the Josephson 
current can be obtained from the formula 19 



/ = -^^tanh( £CT ,„/2fc B T) d£ 



(4) 



where the factor 2 accounts for the valley degeneracy. 

From Eq. ([4]) we have calculated the Josephson criti- 
cal current. Figure [5] shows the dependence of resulting 
Josephson coupling I c Rn on the exchange energy h/Ep 
scaled in units of Fermi energy at zero temperature and 
for EpL/hvF = 10. Here R N = (e 2 /^E„,<T K?)^ 1 is 
the resistance of the corresponding normal (nonsupercon- 
duting) structure. Note that Rn decreases monotonia- 
cally with the exchange field for h > Ep, due to the linear 
increase of the density of states in both spin-subbands. 
In spite of showing the regular cusps form variations indi- 
cating — 7r transitions for all h/Ep, the overall behavior 
of the coupling in the regimes of h < Ep and h > Ep 
is drastically different. For h < Ep the envelope of the 
curve decreases monotonically with h/Ep to reach a min- 
imum for a half-metallic graphene h = Ep, just similar 
to the behavior of I c Rn in a common SFS junction^. 
However for h > Ep the coupling has a finite value and 
surprisingly increases smoothly with h before showing a 
slow damping at strong exchange fields h » Ep. Thus 
I c Rn as a function of h develops a smooth maximum at 
an exchange field which we have found to depend weakly 
on the doping of F graphene EL/hvp. 



We can understand the above behavior of I c Rn in 
terms of the change in the configuration of the two spin 
subbands by varying the ratio h/Ep. For a normally n- 
doped F graphene (Ep > 0), while for h < Ep the charge 
carriers in both spin subbands are of the same n-type 
(Fig. Hk), for h > Ep the spin-down carriers turn into 
the p-type as the Fermi level shifts into the valence sub- 
band (Fig. QJ>). For h < Ep the Andreev bound states 
which carry the supercurrent are made of electrons and 
holes with the same type n as shown in Fig. QJi. The AR 
of the excitations at the Fermi level is of retro type as 
in an ordinary FS interface. In this regime by increasing 
h/Ep the amplitude of AR decreases, which leads to a 
decline in the Josephson coupling of the SFS junction. 
However when h > Ep Andreev correlated electron-hole 
pairs in F are of different n and p types which are cou- 
pled via a specular— ^ AKR at FS interfaces. Due to the 
exchange field induced enhancement of the amplitude of 
AKR, the resulted Josephson coupling I c Rn shows an 
increase by h for h > Ep . But this increase does not con- 
tinue for higher h/Ep where I c Rn decreases very slowly 
after showing a smooth maximum. This slow decline of 
I c Rn, in spite of the increase in the amplitude of AKR, is 
the result of superposition of the contributions of differ- 
ent transverse modes in the supercurrent (see Eq. (j4|)). 

The behavior of Josephson junction with a half metal 
graphene h = Ep in which the density of states of the 
down-spin subband goes to zero is even more dramatic. 
In spite of the fact that there is no propagating excitation 
at the Fermi level with down spin, we still find a nonzero 
critical current for h = Ef as it is seen in Fig. [2l Indeed, 
in this case we have only Andreev bound states combined 
from propagating spin up electrons (holes) and evanes- 
cent spin down holes (electrons). Such Andreev states 
exist for all values of h/Ep, however their contribution 
to the supercurrent is negligible unless at the vicinity of 
h = Ep where they play the main role. 

Now let us analyze the effect of a finite temperature. 
From Eq. (0| we have found that the ballistic graphene 
SFS junction can transit from 0-state to 7r-state by vary- 
ing the temperature. The resulting phase diagram in the 
plane of T/T c and h/Ep is shown in Fig. QJ around 
h = Ep and when EL/hvF = 10. As it can be seen the 
values of exchange fields in which — tt transitions oc- 
cur, increase (decrease) with temperature for the regime 
of h < Ep (h > Ep) such that the phase boundaries for 
h > E are the mirror form of those for h < E. This 
— 7r phase diagram is different from that of a common 
ballistic SFS junction^. The temperature-induced — n 
transition can also be seen from the inset of Fig. [2] where 
we have compared the oscillations of I c Rn at two tem- 
peratures T = and T = 0.99T C . We see that the place 
of a — 7r cusp depends on T. This dependence is differ- 
ent in two regimes of h < Ep and h > E, which results 
in different forms of the phases boundaries as described 
above. Note that as the result of a strongly nonsinu- 
soidal current-phase relation at T = 0, there is a large 
residual supercurrent at a transition point 17 . As T — > T c 
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the current-phase relation becomes pure sinusoidal and 
the residual supercurrents vanish. 

We note that the long range Josephson coupling in 
graphene SFS junctions is carried by the superconduct- 
ing correlations of two electrons with opposite spins. This 
effect, arising from the Dirac-like spectrum of excitations 
and their chiral nature, is fundamentally distinct from 
the recently discovered long range proximity effect in or- 
dinary SFS structures, which was attributed to the gen- 
eration of the spin-parallel triplet correlations by inho- 
mogeneity in the direction of the exchange field2i. 

The practical importance of the effects predicted here 
is connected with the possibility of fabricating high qual- 
ity FS structures in graphene which seems to be quite 
feasible by considering the recent experimental realiza- 
tions of proximity induced superconductivity^^ and fer- 
romagnetic correlation o 23 ' 24 ' 25 in graphene. In addition 
to the proximity induced correlations 2 ^, intrinsic fer- 
romagnetism were also predicted to exist in graphene 
sheets^ and nanoribbons 2 ^. One alternative way to pro- 



duce Josephson F contact would be doping of the spac- 
ing part between two S regions (on top of which metallic 
superconducting electrodes are deposited) by magnetic 
atom impurities 2 ^. 

In conclusion, we have demonstrated the existence of 
a long range supercurrent in weakly doped graphene fer- 
romagnetic Josephson junctions. In contrast to the com- 
mon view, a half metallic graphene shows a nonvanishing 
Josephson coupling I c Rn which increases by increasing 
the exchange field h above the Fermi energy Ep , and 
shows only a slow damping at strong exchanges /i> Ep. 
We have explained this long range coupling as the re- 
sult of the exchange field mediated Andreev-Klein pro- 
cess at FS interfaces, which enhances the induction of 
superconducting correlations between electrons with op- 
posite spins in F. We have also presented the — ir phase 
diagram of the coupling in the plane of T/T c and h/ Ep, 
which reveals the distinct shapes of the phase boundaries 
in two cases of h > Ep and h < Ep. 
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